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We present a study of the process e'*'e“ —>■ using data samples of 1092 pb~^ at 

^/s = 4.23 GeV and 826 pb“^ at = 4.26 GeV collected with the BESIII detector at the 
BEPGII storage ring. With full reconstruction of the D meson pair and the bachelor tt^ in the 
final state, we confirm the existence of the charged structure Zc(3885)^ in the (DD*)^ system in 
the two isospin processes e^e“ —>■ and e^e“ D~. By performing a simul¬ 

taneous fit, the statistical significance of .Zc(3885)^ signal is determined to be greater than lOcr, 
and its pole mass and width are measured to be Mpoie=(3881.7±1.6(stat.)±1.6(syst.)) MeV/c^ and 
rpoie=(26.6±2.0(stat.)±2.1(syst.)) MeV, respectively. The Born cross section times the (DD*)^ 
branching fraction —>■ 7r^Zc(3885)^) x Br{Zc{3885)^ —>■ (DD*)^)) is measured to be 

(141.6 ± 7.9(stat.) ± 12.3(syst.)) pb at = 4.23 GeV and (108.4 ± 6.9(stat.) ± 8.8(syst.)) pb at 
^/s = 4.26 GeV. The polar angular distribution of the 7r^-Zc(3885)^ system is consistent with the 
expectation of a quantum number assignment of = 1^ for ^c(3885)^. 

PACS numbers: 14.40.Pq, 13.25.Gv, 12.38.Qk 


I. INTRODUCTION 

The V(4260) was first observed by BaBar in 
the initial-state-radiation (ISR) process e“''e“ —>■ 
T/Satt+tt" J/'0 Q. This observation was subsequently 
confirmed by CLEO 0 and Belle Q. Unlike other char- 
monium states, such as "0(4040), "0(4160) and 0(4415), 
Y (4260) does not have a natural place within the quark 
model of charmonium [^. Many theoretical interpreta¬ 
tions have been proposed to understand the underlying 
structure of F(4260) [l-Q, more precise experiments are 
necessary to give a decisive conclusion. 

In recent years, a common pattern has been observed 
for the charmoniumlike states in the systems tt J/0, Trip', 
Trhc and nxc as well as in pairs of charmed mesons DD* 
and D*D*. Belle observed some charged structures called 
Z(4430)^ in the system [8l-[Toj|. and Zi(4050)^ and 
^ 2 ( 4250 )^ in the tt^Xci invariant mass spectra [ll| in 
B meson decays. The Z(4430)^ has recently been con¬ 
firmed by LHCb [l^l in the 7r=*=0' system. However, nei¬ 
ther Zi(4050)^ nor Z2(4250)^ are found to be significant 
in BaBar data [Hill. BESIII [3 and Belle 0 ob¬ 
served the Zc(3900)^ in the J/0 invariant mass distri¬ 
bution in a study of e“'"e“ —>■ tt~ J/ip] this observation 

was confirmed with CLEOc data at -0i=4.17 GeV [la. 
More recently, BESIII has reported the observations of 
the Zc(3900)° in the n^J/ip system [I^, Zc(4020) in 
the Trhr system m Z,(4025) in the D*D* 
tem [iJlia, and Zc(3885)^ in the {DD*)^ system Q. 
It is interesting to note that all these states lie close to the 
threshold of some charm meson pair systems and some 
of them even have overlapping widths. It is therefore 
important to obtain more experimental information to 
improve the understanding of all these states. 

In a previous paper by BESIII [l^, a structure 
called Zc(3885)^ was observed in the study of e+e“ —>■ 
{D^ —> K~Tr~^) and e+e“ —>■ 

{D~ K^TT~Tr~) using a 525 pb“^ subset of the data 
sample collected around i/s = 4.26 GeV. That study em¬ 
ploys a partial reconstruction technique by reconstruct¬ 
ing one final-state D meson and the bachelor tt coming 
directly from e+e“ decay (“single D tag”or ST) and in¬ 
ferring the presence of the D* from energy-momentum 


conservation. In this analysis, we present a combined 
study of the processes e+e“ —>■ (tt+D^D^- 

tagged) and e+e“ — D~-tagged) using 
data samples of 1092 pb“^ at y^= 4.23 GeV and 826 pb~^ 
at -05=4.26 GeV [ 2 ^ collected with the BESIII detector 
at the BEPGII storage ring (charge conjugated processes 
are included throughout this paper). We reconstruct the 
bachelor 7r+ and the D meson pair (“double D tag”or 
DT) in the final state. Because the tt from D*~ and 
D*^ decays has low momentum, it is difficult to recon¬ 
struct directly. We denote it as the “missing tt” and infer 
its presence using energy-momentum conservation. The 
£>0 mesons are reconstructed in four decay modes and 
the D~ mesons in six decay modes. The double D tag 
technique allows the use of more D decay modes and 
effectively suppresses backgrounds. 


II. EXPERIMENT AND DATA SAMPLE 

The BESIII detector is described in detail else¬ 
where [ 2 ^. It has an effective geometrical acceptance 
of 93% of 47r. It consists of a small-cell, helium-based 
(40% He, 60% CaHg) main drift chamber (MDC), a plas¬ 
tic scintillator time-of-flight system (TOE), a CsI(TI) 
electromagnetic calorimeter (EMC) and a muon system 
(MUG) containing resistive plate chambers (RPC) in 
the iron return yoke of a 1 T superconducting solenoid. 
The momentum resolution for charged tracks is 0.5% at 
a momentum of 1 GeV/c. Charged particle identifica¬ 
tion (PID) is accomplished by combining the energy loss 
(dE/dx) measurements in the MDC and flight times in 
the TOP. The photon energy resolution at 1 GeV is 2.5% 
in the barrel and 5% in the end caps. 

The GEANT4-based 1^. 1^ Monte Carlo (MC) sim¬ 
ulation software BOOST includes the geometric and 
material description of the BESIII detectors, the detec¬ 
tor response and digitization models, as well as the track¬ 
ing of the detector running conditions and performance. 
It is used to optimize the selection criteria, to evaluate 
the signal efficiency and mass resolution, and to estimate 
the physics backgrounds. The physics backgrounds are 
studied using a generic MC sample which consists of the 
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production of the F(4260) state and its exclusive de¬ 
cays, the process e+e” —>■ D^*\ the production 

of ISR photons to low mass ip states, and QED pro¬ 
cesses. The y(4260) resonance, ISR production of the 
vector charmonium states, and QED events are gener¬ 
ated by KKMC [1^. The known decay modes are gen¬ 
erated by EVTGEN [Ullll with branching ratios being 
set to world average values from the Particle Data Group 
(PDG) [^ . and the remaining unknown decay modes are 
generated by LUNDCHARM In addition, exclu¬ 

sive MG samples for the process e+e“ — DjD*,Dj —>• 
d(*) 7r(7r) are generated to study the possible background 
contributions from neutral and charged highly excited D 
states (denoted as Dj, where J is the spin of the meson), 
such as D^{2A00), Di(2420), Di(2430) and £>^(2460). To 
estimate the signal efficiency and to optimize the selec¬ 
tion criteria, we generate a signal MG sample for the 
process e+e” 7r+Ze(3885)- (^^(3885)- ^ {DD*)-) 
and a phase space MG sample (PHSP MG) for the pro¬ 
cess e+e" —?► 7r+(DZ)*)“. Here the spin and parity of the 
Zc(3885)“ state are assumed to be 1+, which is consis¬ 
tent with our observation. 


III. EVENT SELECTION AND BACKGROUND 
ANALYSIS 

Charged tracks are reconstructed in the MDC. For 
each good charged track, the polar angle must satisfy 
I cos0| < 0.93, and its point of closest approach to the 
interaction point must be within 10 cm in the beam 
direction and within 1 cm in the plane perpendicular 
to the beam direction. To assign a particle hypothesis 
to the charged track, dE/dx and TOE information are 
combined to form a probability Prob(Ar) (Prob(7r)). A 
track is identihed as a AT (tt) when Prob(Ar) > Prob(7r) 
(Prob(7r) > Prob(Ar)). Tracks used in reconstructing Kg 
decays are exempted from these requirements. 

Photon candidates are reconstructed by clustering 
EMC crystal energies. For each photon candidate, the 
energy deposit in the EMC barrel region (| cos0| < 0.8) is 
required to be greater than 25 MeV and in the EMC end- 
cap region (0.84 < | cos0| < 0.92) greater than 50 MeV. 
To eliminate showers from charged particles, the angle 
between the photon and the nearest charged track is re¬ 
quired to be greater than 20°. Timing requirements are 
used to suppress electronic noise and energy deposits in 
the EMC unrelated to the event. 

We reconstruct 7r° candidates from pairs of photons 
with an invariant mass in the range 0.115 < M^-y < 
0.150 MeV/c^. A one-constraint (1C) kinematic fit is 
performed to improve the energy resolution, with 
being constrained to the known mass from PDG [32j| . 

Kg candidates are reconstructed from pairs of oppo¬ 
sitely charged tracks which satisfy |cos0| < 0.93 for 
the polar angle and the distance of the track to the 
interaction point in the beam direction within 20 cm. 


For each candidate, we perform a vertex fit constraining 
the charged tracks to a common decay vertex and use 
the corrected track parameters to calculate the invari¬ 
ant mass which must be in the range 0.487 < M^+^- < 
0.511 GeV/c^. To reject random tt+tt” combinations, a 
secondary-vertex fitting algorithm is employed to impose 
a kinematic constraint between the production and decay 
vertices [sj. 

The selected 7r=*=, K^, 7r° and Kg are used to recon¬ 
struct D meson candidates for the and D~D^ 

double tag. The candidates are reconstructed in 
four final states: K~tt~^, K~tt~^tt^, and 

(in the following labeled as 0, 1, 2, 
and 3, respectively), and the D~ candidates in six fi¬ 
nal states: AT+tt^tt”, K'^ tt~'k~ ^ KgTT~, AT^Tr^Tr^, 
KgTT^TT~TT~ and K^K~'k~ (labeled as A, B, C, D, E 
and E, respectively). If there is more than one candi¬ 
date per possible DT mode, the candidate with the min¬ 
imum AM is chosen, where AM is the difference be¬ 
tween the average mass M = [M{D) + M{D)]/2 and 
[Mpdg(.D)-|-Mpdg(.D)]/2 (Mpdg(.D) andMpDG(.D) are 
the D mass and D mass from PDG [s^, respectively). 
Figure [T] shows the distributions of M{D) versus M{D) 
for all DT candidates at -^5=4.26 GeV. The combina¬ 
torial background tends to have structure in AM but is 
flat in the mass difference AM = M{D) — M{D). The 
signal region in the M{D) versus M{D) plane is defined 
as -20 < AM < 15 MeV/c^ (-17 < AM < 14 MeV/c^) 
and |AM| < 40 MeV/c^ (|AM| < 35 MeV/c^) for D°D° 
{D~D°) candidates. 

To reconstruct the bachelor 7r+, at least one addi¬ 
tional good charged track which is not among the de¬ 
cay products of the D candidates is required. To re¬ 
duce background and improve the mass resolution, we 
perform a four-constraint (4C) kinematic fit to the se¬ 
lected events. It imposes momentum and energy conser¬ 
vation, constrains the invariant mass of D (D) candidates 
to Mpdg(A>) {Mpug{D)), and constrains the invariant 
mass formed from the missing tt and the corresponding 
D candidate to Mpdg(.D*) [^. This gives a total of 
7 constraints. The missing tt three-momentum needs to 
be determined, so we are left with a four-constraint fit. 
The of the 4C kinematic fit (x^c) required to be 
less than 100. If there are multiple candidates in an 
event, we choose the one with minimum Xac- sup¬ 
press the background process e+e" —>■ D*D*, we require 
M{-k+D^)> 2.03 GeV/c2 (M(7r+D-) > 2.08 GeV/c^) 
for 7r+£>°£>°-tagged (7r+£)“A)°-tagged) events. We de¬ 
fine the reconstructed Dtt recoil mass Mj-ecoiiiDn) via 

Mrecoil(A’7r)2c4 = (^Ecm- Ed - E^)'^ - \pcm-PO, 

where (Acm, Pcm), [Ed, Pd) and [E^,, p^) are the four 
momentum of the e^e“ system, D and tt in the e“'"e“ 
rest frame, respectively. Figure [5] shows the Mrecoii(A7r) 
distributions at ^/s=4.26 GeV after all of the above selec¬ 
tion criteria. The results of signal MG and PHSP MG are 
provided to verify the signal processes and optimize the 
selection criteria. A study of generic MC sample shows 
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FIG. 1. Masses of the D and D candidates for all DT modes at y^=4.26 GeV. The vertical (horizontal) bands centered at 
M{D) {M{D)) contain the DT candidates in which the D (D) candidate was reconstructed correctly, but the D (D) was not. 
The diagonal bands contain the “mis-reconstructed” DD candidates (all of the D and D final states were reconstructed, but one 
or more final states from the D were interchanged with corresponding particles from the D). Other combinatorial candidates 
with minimum AM also spread along the diagonal. The left plot shows M{D^) versus while the right plot shows 

M{D~) versus M{D^). The solid rectangles show the signal regions. 


that very few background events which can satisfy the 
above requirements. 

To select the ttDD* events, we require that 
|.Wrecoii(T^7r) — Mpj:,g{D*)\ < 30 MeV/c^. After im¬ 
posing all of the above requirements, a peak around 
3890 MeV /(? is clearly visible in the kinematically con¬ 
strained DD* mass distributions for selected 

events, as shown in Fig.[3l For the 7r+II“I?°-tagged pro¬ 
cess, some events from the isospin partner decay channel 
e+e” —>■ {D*~ —?► can satisfy the above 

requirements, but with different reconstruction efficiency 
and mass resolution. We treat these as signal events and 
combine them with the 7r+II”£)°-tagged process. For 
the data sample at -^5=4.23 GeV, we employ the same 
event selection criteria and obtain similar results. 

We use the generic MC sample to investigate pos¬ 
sible backgrounds. There is no similar peak found 
near 3.9 GeV/c^ and the selected events predominantly 
have the same final states as Tr~^{DD*)~. From a 
study of the Monte Carlo samples of highly excited 
D states, we conclude that only the process e+e" —>■ 
IIi(2420)Z), I?i(2420) —^ ttD* can produce a peak near 
the threshold in the DD* mass distribution, although 
the probability of this is small due to the kinematic 
boundary. To examine this possibility, the events are 
separated into two samples according to | cosSttuI <0.5 
and Icos^Trul > 0.5, where is the angle between 
the directions of the bachelor 7r+ and the D meson in 


the DD* rest frame. Defining the asymmetry A = 
in>o.5 - »^<o.5)/(?t->o. 5 + »^<o.5), where n>o.5 and n<o.5 
are the numbers of events in each sample, we found that 
the asymmetry in data, Mdato=0.11±0.07, is compati¬ 
ble with the asymmetry expected in signal MC, A^f^ 
=0.01±0.01, and incompatible with the expectations for 
DDi(2420) MC, =0.43±0.01. Considering the 

kinematic boundary of this process, we conclude that the 
ZIZ)i(2420) contribution to our observed Born cross sec¬ 
tion is smaller than its relative systematic uncertainty. 
This is consistent with the ST analysis [^ . 


IV. SIGNAL EXTRACTION 


To extract the resonance parameters and yield of 
Zc(3885)“ in the {DD*)~ mass spectrum, both processes 
are fitted simultaneously with an unbinned maximum 
likelihood method using two different data samples at 
Y^=4.23 GeV and i/i=4.26 GeV. The {DD*)~ invariant 
mass distribution is described as the sum of two proba¬ 
bility density functions (PDFs) representing the signal 
and background. The signal PDF is given by 


PDF(m£,£,.) 


0 A]e(TOji^.) 


( 1 ) 


where the integral is performed over the fit range of the 
{DD*)~ mass spectrum, S'(m£i£i.) is the signal term 
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FIG. 2. The Mj-ecoii distributions for (a) 7r'''D°_D°-tagged events and (b) 7r'''D“II°-tagged events at y^=4.26 GeV. The 
dots with error bars are data. The dashed (red) and solid (blue) lines are signal MC and PHSP MC, respectively. The arrows 
(pink) indicate nominal selection criteria. 




FIG. 3. Simultaneous fits to the M{DD*) distributions of ((a) and (c)) 7r^D°I)°-tagged and ((b) and (d)) 7r’''D~D°-tagged 
processes for ((a) and (b)) data at y^=4.23 GeV and for ((c) and (d)) data at y^=4.26 GeV. The dots with error bars are 
data and the lines show the projection of the simultaneous fit to the data. The solid lines (blue) describe the total fits, the 
dashed lines (red) describe the signal shapes and the green areas describe the background shapes. 
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convolved with the mass resolution, and ) is the 

reconstruction efficiency. The background PDF is pa¬ 
rameterized by phase space MC simulation. The sig¬ 
nal and background yields and the mass and width of 
Zc(3885)“ are determined in the fit. The mass and width 
of Zc(3885)“ are constrained to be the same for both pro¬ 
cesses. 


A. Signal Term 

The process e+e“ —>■ 7r+Zc(3885)“ with Zc(3885)“ —^ 
I is described with phase space generalized for the angular 
momentum L of the 7r+ — Zc(3885)“ system, where I 
denotes D~D*^ (labeled as a) and D^D*~ (labeled as 
h). The Zc(3885)“ is described by a mass dependent 
width Breit-Wigner (MDBW) parameterization 

Si{moD*) oc dN/dmoD* 


where k* is the momentum of Zc(3885) in the e“'"e rest 
frame, /l(k*) is the Blatt-Weisskopf barrier factor [s^ . 


oc 




TTb) 


,( 3 ) 


Fi = TzM/<li?^~^^[mzJrnjjQ,][fe{q^)/fe{q^)]‘^, gf is 
the D momentum in the Zc(3885)“ rest frame, i is the 
angular momentum of the {DD*)~ system, and = 
ql{mz^). In the fit, mz^ and Tz^ are free parameters, 
while L = 0 and £ = 0 are fixed according to the analysis 
of angular distributions below. Parameters of the reso¬ 
lution and efficiency functions, obtained from MC and 
described below, are fixed in the fit. 


B. Reconstruction Efficiency and Mass Resolution 

In order to obtain the reconstruction efficiency and 
mass resolution, we generate a set of MC samples for 
e+e“ —>• {Z~ —>• {DD*)~), each with a fixed mass 

value, zero width and = I’*' of the Z~, and subject 
these MC samples to the same event selection criteria. 
The isospin channel e+e” —>■ {D*~ —>■ D~Tr^) 

can feed into the 7r+D“D°-tagged process. We there¬ 
fore generate two corresponding MC samples by assum¬ 
ing the same decay branching fraction between the pro¬ 
cess Z~ —>• D~D*° and Z~ —^ D°D*~. The reconstruc¬ 
tion efficiency is estimated using the sum of the two MC 
samples, as shown in Fig. |4l 

MC samples for e+e“ —>■ tt^Z~ {Z~ —>• {DD*)~) are 
used to determine the mass resolution. The mass and 
width of Zc are set to be 3890 MeV/c^ and 0 MeV, re¬ 
spectively. The mass resolution for the 7r~*~T)°Z)°-ta gged 
process is described by a Crystal Ball (CB) function [3^ . 
Since the 7r+D“D°-tagged process contains two isospin 
processes, the mass resolution is represented by a sum 


of two CB functions with a common mean and differ¬ 
ent widths. The fit results for both processes are shown 
in Fig. [5l The resolution for the 7r+Z)°Z)°-tagged pro¬ 
cess is determined by the fit to be I.liO.l MeV/c^, 
while the resolution for the 7r+D“D°-tagged process is 
calculated to be 2.2±0.1 MeV/c^ using the equation 
/iCTi + (1 — where cti and CT 2 are the individual 

widths of each of the two CB functions and fi is the 
fractional area of the first CB function. 


C. Fit Results 

As shown in Fig. [3l we perform a simultaneous fit 
to the M{DD*) distributions for the 7r+D°Z)°-tagged 
and 7r+D”D°-tagged processes with •y/s=4.23 GeV and 
i/i=4.26 GeV data samples. The statistical signifi¬ 
cance of Z'c(3885)“, estimated by the difference of log- 
likelihood values with and without signal terms in the fit, 
is greater than lOcr. The mass and width of Zc(3885)“ 
are fitted to be Mz^[ 3 S 85 ) = (3890.3±0.8) MeV/c^ and 
Fz,,( 3885 ) = (31.5±3.3) MeV, where the errors are statis¬ 
tical only. Since the resulting mass and width might 
be different from the actual resonance properties due 
to the parameterization function of .Zc(3885), we cal¬ 
culate the pole position (P = Mpoie — irpoie/2) of 
.Zc(3885) which is the complex number where the de¬ 
nominator of BWi{muQ,) is zero, and regard Mpoie 
and Fpoie as the final result. The corresponding pole 
mass (Mpoie) and width (Fpoie) of Zc(3885) are Mpoie 
= (3881.7±1.6) MeV/c^ and Fpoie = (26.6±2.0) MeV, 
respectively. 


D. Angular Distribution 

The quantum number assignment for Zc(3885)“ 
is investigated by examining the distribution of |cos0,r|, 
where 0,r is the tt’*' polar angle relative to the beam di¬ 
rection in the center-of-mass frame. If = I'*', the 
relative orbital angular momentum of the 7r“'"-Zc(3885)“ 
system could be either S'-wave or D-wave. If we neglect 
the small contribution of D-wave due to the closeness of 
the threshold, the |cos0t| distribution is expected to be 
flat. If = 0“ (1“), the 7r+-Zc(3885)“ system occurs 
via a P-wave and the |cos0,r| is expected to follow sin^0^ 
(l-|-cos^0,r) distribution. 

The |cos0^| distribution of data is plotted with the ef¬ 
ficiency corrected signal yield of combined data samples 
at -^5=4.23 GeV and -^5=4.26 GeV in ten |cos 07 r| bins, 
where the signal yields in different bin are extracted 
with the same simultaneous fit method described above. 
Figures [ 5 ] (a) and (b) show the |cos0,r| distribution for 
7r+D°D°-tagged process and 7r+D“ZI°-tagged process, 
respectively. The data agrees well with the flat distribu¬ 
tion expected for = 1+ (y^/NDF = 16.5/9 for the 
7r+II°i)°-tagged process and 12.8/9 for the 7r~^D~D^- 
tagged process) and disagrees with the sin^0^ distribu- 
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FIG. 4. Distributions of the efficiency versus M{DD*) for ((a) and (c)) 7r"''D°l)°-tagged and ((b) and (d)) 7r"*'D~D°-tagged 
processes at ((a) and (b)) y^=4.23 GeV and ((c) and (d)) ■\/s=4.26 GeV. The dots with error bars are the efficiencies determined 
from MG. The curves show the fits with a piecewise linear function. 


tion expected for = 0“ (x^/NDF = 103.1/9 for the 
7 r+£)°Z)°-tagged process and 104.9/9 for the 'k^D~D^- 
tagged process) and = 1~ (x^/NDF = 106.3/9 for the 
7 r+£)°Il°-tagged process and 104.9/9 for the 
tagged process), where NDF is the number of degrees of 
freedom in the fit. 


E. Born Cross Section 

For the 7r+Il°Z)°-tagged process, the Born cross sec¬ 
tion times the {DD*)~ branching fraction of Zc(3885)“ 

(cr X Br) can be calculated by 

cr(e+e- ^ 7r±Zc(3885)=F) x Br(Zc(3885)=F [DD*)^) 

^ _ N __( 4 ) 

£(1 -I- i5’’)(l -I- S'") j ei^BriBr^BriD*- ir-P^)! \ 


where N is the signal yield, C is the integrated lumi¬ 
nosity, eij is the signal efficiency for the 7r+Il°Z)°-tagged 
process listed in Table Hill of Appendix where the sub¬ 
scripts i,j = 0...3 denote the neutral D final state, 
Bvi is the individual branching fraction for D decay 
from PDG [1^, the radiative correction factor (1 -|- S'") 
is determined by the measurement of the line shape of 
cr(e“''e“ —>■ ttDD*) [2^, the vacuum polarization factor 
(1 -I- (i") is considered in the MC simulation [s^ and / = 
Br{Zc{3885)- D°D*-)/Br{Zc{3885)- (DD*)-) 

= 0.5, assuming isospin symmetry. The value of all above 
variables are listed in Table HI 

Since the 7r+T)”Il°-tagged process contains two pro¬ 
cesses of Zc{3885)- D-D*° with D*° tt°D° (la¬ 

beled as a) and 2'c(3885)“ — ?► D^D*~ with D*~ —>■ Tr^D~ 
(labeled as /3), the Born cross section times the {DD*)~ 
branching fraction of Zc(3885)“ can be given by 


cr(e+e“ ^ 7r±Zc(3885)=F) x Br(Zc(3885)=F ^ (DD*)^) 

N 

£(1 -I- (5’')(1 -I- e'^jBriBrjBr{D*° -)> Tr°D°) + Y,ij eijBriBrjBr{D*- 


( 5 ) 
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FIG. 5. Fits to the mass resolution at 3890 MeV for ((a) and (c)) 7r'''D°I)°-tagged and ((b) and (d)) 7r"*'H”D°-tagged processes 
at ((a) and (b)) v^=4.23 GeV and ((c) and (d)) v^=4.26 GeV. The dots with error bars show the distributions of mass 
resolutions obtained from MC, the curves show the fits. 


where eg and e^- are the signal efficiency for the two 
TT^D~-tagged processes, listed in Table IIVI and m 
of Appendix El the subscripts i and j denote the D 
and final states, respectively, with i = A... F and 
j = 0...3, Br{D*° Tr°D°) = (61.9±2.9)% and 
Br{D*- tt°D-) = (30.7±0.5)% [3|. The value of 
all above variables are listed in Table HI 


We also add a Zc(4020) in the fit with mass and width 
fixed to the BESIII measurement M- The fit prefers the 
presence of a Zc(4020)“ with a statistical significance of 
l.Ocr. We determine the upper limit on cr x Br at the 
90% confidence level (C.L.), where the probability den¬ 
sity function from the fit is smeared by a Gaussian func¬ 
tion with a standard deviation of the relative systematic 
error in the a x Br measurement. We obtain a{e'^e~ —>■ 
7 r±Zc( 4020 )=F) X Br(Ze(4020)=F ^ (DD*)^) < 18 pb at 
•^5=4.23 GeV and <15 pb at •v/i=4.26 GeV, respectively, 
at 90% C.L.. 


V. SYSTEMATIC UNCERTAINTIES 

The systematic uncertainties for the pole mass and 
width of Zc( 3885)“, and the product of Born cross sec¬ 
tion times the {DD*)~ branching fraction of Zc(3885)“ 
(ct X Br) are described below and summarized in Table HIl 
The total systematic uncertainty is obtained by summing 
all individual contributions in quadrature. 

Beam Energy: In order to obtain the systematic uncer¬ 
tainty related to the beam energy, we repeat the whole 
analysis by varying the beam energy with ± 1 MeV in the 
kinematic fit. The largest difference on the pole mass, 
width and the signal yields is taken as a systematic un¬ 
certainty. 

Mass Calibration: The uncertainty from the mass cali¬ 
bration is estimated with the difference between the mea¬ 
sured and nominal D* masses. We fit the D* mass spec¬ 
tra calculated with the output momentum of the kine¬ 
matic fit described in the Sec. uni after removing the D* 
mass constraint. The deviation of the resulting D* mass 
to the nominal values is found to be 0.84±0.16 MeV/c^. 
The systematic uncertainty due to the mass calibration 
is taken to be 1.0 MeV/c^. 
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FIG. 6. Fits to IcosSti-I distributions for (a) 7r’'"D'’l)'’-tagged and (b) 7r^Z)“D°-tagged processes. The dots with error bars show 
the combined data corrected for detection efficiency at y^=4.23 GeV and -^5=4.26 GeV, the solid lines show the fits using 
= 1'*" hypothesis, and the dashed and dotted curves are for the fits with = 0“ and — 1“ hypothesis, respectively. 


TABLE I. Summary of the product of Born cross sections times the {DD*) branching fraction of ^c(3885) (cr x Br), the 
errors are statistical only. 



7r"''D’^D’^-tagged process 

D’^-tagged process 

4.23 GeV 

4.26 GeV 

4.23 GeV 

4.26 GeV 

N 

384±30 

207±18 

418±34 

239±22 

B (pb-^) 

1091.7 

825.7 

1091.7 

825.7 

l+<5" 

0.89 

0.92 

0.89 

0.92 


1.056 

1.054 

1.056 

1.054 

a X Br (pb) 

147.5±11.5 

109.2±9.7 

136.6±11.0 

107.5±9.7 


£(l + 5'’)(l + 5’') : The integrated luminosities of the 
data samples are measured using large angle Bhabha 
events, with an estimated uncertainty of 1.0% [2^. The 
systematic uncertainty of the radiative correction factor 
is estimated by changing the parameters of the line shape 
of CT(e+e“ —>■ ttDD*) within errors. We assign 4.6% as 
the systematic uncertainty due to the radiative correction 
factor according to Ref. . The systematic uncertainty 
of the vacuum polarization factor is 0.5% [s^. 

Signal shape: The systematic uncertainty associated 
with the Zc(3885)“ signal shape is evaluated by repeating 
the fit on the M{DD*) distribution with a mass constant 
width BW line shape (MCBW, —^- 2 —-—^ r—) for 

Zc( 3885)“ signal. The resulting difference to the nominal 
results are taken as a systematic uncertainty. 

Zc(4020)“ signal The systematic uncertainty associated 
with the possible existence of the Zc(4020)“ in our data 
is estimated by adding the Z’c(4020)“ in the fit. The dif¬ 
ference of fit results is taken as a systematic uncertainty. 

Background shape: The systematic uncertainty due to 
the background shape is investigated by repeating the fit 
with function fbkgimaB,) oc - M^znYiM^ax - 

'^dd*Y 0 foi' the background line shape, where 
and Mmax are the minimum and maximum kinematically 


allowed masses, respectively, c and d are free parameters. 
The resulting difference to the nominal results is taken 
as a systematic uncertainty. 

Fit bias: To assess a possible bias due to the fitting 
procedure, we generate 200 fully reconstructed data-size 
samples with the parameters set to the values (input val¬ 
ues) returned by the fit to data. Then we fit these sam¬ 
ples using the same procedures as we fit the data, and 
the resulting distribution of every fitted parameter with 
a Gaussian function. The difference between the mean 
value of Gaussian and the input value is taken as a sys¬ 
tematic uncertainty of the ht bias. 

Signal region of DT: In order to obtain the systematic 
uncertainty related to the selection of the signal region 
of double D tag, we repeat the whole analysis by chang¬ 
ing the signal region in the M{D) versus M{D) plane 
from the nominal region to —15 < AM < 10 MeV/c^ 
(|AM| < 30 MeV/c^) and -25 _< AM < 20 MeV/c^ 
(lAMj < 60 MeV/c^) for 7r+Il°Il°-tagged, and —14 < 
AM < 11 MeV/c^ (|AM| < 28 MeV/c^) and -20 < 
AM < 17 MeV/c^ (|AM| < 42 MeV/c^) for tt+D-D°- 
tagged processes. The largest difference of fit results is 
taken as a systematic uncertainty. 

Efficiency related: We refer to the system- 
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atic uncertainty for j ^ijBriBrjBr{D*~ —>• 

Tr~D°) and efjBriBrjBr{D*^ —>• 7r°_D°) + 

j ^ijBriBrjBr{D*~ —>• tt°D~)) as the efficiency 
related systematic uncertainty for 7r+Zl°Z)°-tagged and 
7r+£>“£>°-tagged processes, respectively. The efficiency 
related systematic uncertainty includes the uncertain¬ 
ties from MC statistics, PID, tracking, 7r° and Kg 
reconstruction, kinematic fit, cross feed and branching 
fractions of D and D* decay. The uncertainty due to 
finite MC statistics is taken as the uncertainty of the sig¬ 
nal efficiency. A systematic uncertainty of 1% is assigned 
to each track for the difference between data and sim¬ 
ulation in tracking or PID [2^. For 7r° reconstruction, 
the corresponding uncertainty is 3% per 7r° [s^. For 
Kg reconstruction, the corresponding uncertainty is 4% 
per Kg [ 4 ^. The uncertainty due to the kinematic fit is 
estimated by applying the track-parameter corrections 
to the track helix parameters and the corresponding 
covariance matrix for all charged tracks to obtain im¬ 
proved agreement between data and MC simulation [4lj . 
The difference between the obtained efficiencies with 
and without this correction is taken as the systematic 
uncertainty for the kinematic fit. The cross feed among 
different decay modes is estimated using the signal MC 
simulation as detailed in Tables IVIHVIIII of Appendix |B] 
The systematic uncertainties for the branching fractions 
of D and D* decay are estimated by PDG [13. A 
summary of the systematic uncertainties for signal 
efficiency is listed in Tables IVUlVIIIl of Appendix |B] 
The total efficiency related systematic uncertainties are 
combined by considering the correlation of uncertainties 
between each decay channels. 


VI. SUMMARY 

In summary, based on the data samples of 1092 pb“^ 
taken at -^5=4.23 GeV and 826 pb“^ taken at 
■^5=4.26 GeV, we perform a study of the process 
e+e" —>■ 7 r~ {DD*)~^ and confirm the existence of the 
charged charmoniumlike state Zc(3885)“ in the {DD*)~ 
system. The angular distribution of the 7r+ — Zc(3885)“ 
system is consistent with the expectation from a 
= 1+ quantum number assignment. We perform 
a simultaneous fit to the {DD*)~ mass spectra for 
the two isospin processes of e+e" —>■ Tr'^D^D*~ and 
e+e" —>■ using a mass-dependent Breit 

Wigner function. The statistical significance of the 
Zc(3885) signal is greater than lOcr. The pole mass 
and pole width of Zc(3885)“ are determined to be 
Afpoie=(3881.7±1.6(stat.)±1.6(syst.)) MeV/c^ and 
rpoie=(26.6±2.0(stat.)±2.1(syst.)) MeV, respectively. 
The products of Born cross section and the DD* 
branching fraction of Zc(3885)“ for e“'"e“ —>■ Tr~^D°D*~ 
and e+e” —>■ Tr'^D~D*^ are combined into a weighted 
average [13. For the data samples at -^5=4.23 GeV, the 


result is cr(e“''e —>■ 7r^Zc(3885)^) x Br{Zc{3885)^ —?> 
(DD*)^) = (141.6±7.9(stat.)±12.3(syst.)) pb. 

For the y^= 4.26 GeV data sample, the result is 
CT(e+e- 7r±Zc(3885)=F) x Br{Zc{3885)^ {DD*)^) 
= (108.4±6.9(stat.)±8.8(syst.)) pb. 

The pole mass and pole width of Zc(3885)“ and 
CT(e+e- 7r±Zc(3885)=F) x Br{Zc{3885)^ {DD*)^) 

are consistent with but more precise than those of BE- 
SIIFs previous results [13, with significantly improved 
systematic uncertainties. The improvement in the re¬ 
sults obtained in this analysis is due to the fact that 
the double D tag technique and more D tag modes 
are used and two isospin processes e+e“ —>■ 
are fitted simultaneously with datasets at y/s = 4.23 
and 4.26 GeV. This analysis only has ^9% events in 
common with the ST analysis [13, so the two anal¬ 
yses are almost statistically independent and can be 
combined into a weighted average [13. The combined 
pole mass and width are Mpoie = (3882.2 ± I.I(stat.) ± 
1.5(syst.)) MeV/c^ and Fpoie = (26.5 ± I.7(stat.) ± 
2.1(syst.)) MeV, respectively. The combined CT(e+e“ —>■ 
7r±Zc(3885)=F) x Br{Zc{3885)^ (DD*)^) is (104.4 ± 

4.8(stat.) ± 8.4(syst.)) pb at •y/s=4.26 GeV. 
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TABLE II. Summary of systematic uncertainties on the pole mass and pole width of the 2^(3885)“, and the product of Born 
cross section times the {DD*)~ branching fraction of Zc(3885)~ (cr x Br). The items noted with * are common uncertainties, 
and other items are independent uncertainties. 


„ AMpole ATpole -n-n-axBr V - - 

Source tt+D D -tagged process n^D D -tagged process 


(MeV/c-^) 

(MeV) 4.23 GeV 

4.26 GeV 

4.23 GeV 

4.26 GeV 

Beam Energy 

1.0 

1.6 

3.3 

3.0 

4.9 

3.4 

Mass calibration 

1.0 






£(l-kJ")(l-f 5")* 



4.7 

4.7 

4.7 

4.7 

Signal shape 

0.1 

0.1 

0.1 

0.1 

0.1 

0.1 

Zc(4020)“ Signal 

0.4 

1.0 

2.9 

2.0 

2.8 

3.9 

Background shape 

0.4 

0.1 

2.0 

0.5 

2.9 

0.9 

Fit bias 

0.2 

0.1 

0.5 

0.3 

0.1 

0.8 

Signal region of DT 

0.2 

0.7 

4.2 

1.4 

0.8 

1.4 

Efficiency related 



8.3 

8.3 

7.9 

7.9 

Total 

1.6 

2.1 

11.5 

10.3 

11.2 

10.7 
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Appendix A: Signal Efficiency 

The signal efficiency for 7r+i9°i)°-tagged process at 
•^5=4.23 GeV and i/s=4.26 GeV are listed Table Hill 
while the signal efficiency for 7r+Il“Zl°-tagged process 
and its isospin channel are listed in Table IIVI and |Vl 


Appendix B: The Efficiency Related Systematic 
Uncertainty 

The systematic uncertainties for signal efficiency are 
listed in Table IVIHVIIII 
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TABLE III. Signal efficiency etj (%) for 7r+Zo(3885)“ (Zc(3885)" D°D*- 

denote the neutral D final states: K~'k^, K~'k'^'k^, and K~n 


D* 


TT D°, D'^ —>■ i, D'^ —>■ j, where i and j 
(labeled as 0, 1, 2, 3, respectively). 




=.0 


{*> j} 

0 


1 

2 


3 

4.23 GeV 

4.26 GeV 

4.23 GeV 

4.26 GeV 

4.23 GeV 

4.26 GeV 

4.23 GeV 4.26 GeV 

0 

1 

2 

3 

30.23±0.17 

15.23±0.12 

17.42±0.13 

6.64±0.08 

30.30±0.17 

15.47±0.12 

17.33±0.13 

6.62±0.08 

14.68±0.12 

6.65±0.08 

7.50±0.09 

2.26±0.05 

14.76±0.12 

6.52±0.08 

7.45±0.09 

2.29±0.05 

17.54±0.13 

7.80±0.09 

8.01±0.09 

2.41±0.05 

17.53±0.13 

7.80±0.09 

8.00±0.09 

2.30±0.05 

6.50±0.08 6.46±0.08 
2.45±0.05 2.33±0.05 
2.30±0.05 2.30±0.05 
0.35±0.02 0.30±0.02 


TABLE IV. Signal efficiencies for 7r’'"Zc(3885) (.Zc(3885) —>■ D D*^), —>■ , D i, —>■ j, where i denotes 

the charged D final states: A'’^ 7 r~ 7 r“, A'’^ 7 r~ 7 r“ 7 r°, K%'k~, , K%'k'^'k~'k~ and K'^K~tt~ (labeled as A, B, C, D, E 

and F, respectively), and j denotes the neutral D final states: K~'k^, and A'“ 7 r’'" 7 r"'" 7 r“ 7 r° (labeled as 

0, 1, 2, 3, respectively). 



0 


1 

2 

1 


1 

JS 

4.23 GeV 

4.26 GeV 

4.23 GeV 

4.26 GeV 

4.23 GeV 

4.26 GeV 

4.23 GeV 

4.26 GeV 

A 

24.29±0.16 

23.96±0.15 

11.49±0.11 

11.63±0.11 

13.61±0.12 

13.57±0.12 

4.76±0.07 

4.58±0.07 

B 

10.78±0.10 

10.72±0.10 

4.44±0.07 

4.44±0.07 

4.92±0.07 

4.89±0.07 

1.21±0.03 

1.14±0.03 

C 

24.66±0.16 

25.11±0.16 

12 . 02 ± 0.11 

12.05±0.11 

14.22±0.12 

14.27±0.12 

5.09±0.07 

4.89±0.07 

D 

11.56±0.11 

11.55±0.11 

4.85±0.07 

4.87±0.07 

5.79±0.08 

5.62±0.07 

1.61±0.04 

1.53±0.04 

E 

14.56±0.12 

14.75±0.12 

6.23±0.08 

6.31±0.08 

6.31±0.08 

6.24±0.08 

1.70±0.04 

1.59±0.04 

F 

19.29±0.14 

19.13±0.14 

9.05±0.10 

9.11±0.10 

10.67±0.10 

10.64±0.10 

3.51±0.06 

3.38±0.06 


TABLE V. Signal efficiencies for 7r^Zc(3885) (Zc(3885) —^ D^D* ), D* —>■ tt^D , D i, D° —>■ j, where i and j are 
described in the caption of Table nYi 


U i\ 

0 


1 

2 

1 


1 

Jj 

4.23 GeV 

4.26 GeV 

4.23 GeV 

4.26 GeV 

4.23 GeV 

4.26 GeV 

4.23 GeV 

4.26 GeV 

A 

23.57±0.15 

23.65±0.15 

11.32±0.11 

11.42±0.11 

13.22±0.11 

13.09±0.11 

4.75±0.07 

4.68±0.07 

B 

10.83±0.10 

10.49±0.10 

4.34±0.07 

4.34±0.07 

4.86±0.07 

4.76±0.07 

1.17±0.03 

1.16±0.03 

C 

24.51±0.16 

24.37±0.16 

11.94±0.11 

11.91±0.11 

13.98±0.12 

13.87±0.12 

4.96±0.07 

4.93±0.07 

D 

11.34±0.11 

11.30±0.11 

4.68±0.07 

4.83±0.07 

5.67±0.08 

5.46±0.07 

1.58±0.04 

1.47±0.04 

E 

14.04±0.12 

14.17T0.12 

6.19±0.08 

6.04±0.08 

6.11±0.08 

6.08±0.08 

1.60±0.04 

1.52±0.04 

F 

18.89±0.14 

18.79T0.14 

9.03±0.10 

9.08±0.10 

10.42±0.10 

10.37±0.10 

3.35±0.06 

3.44±0.06 


TABLE VI. The systematic uncertainties for signal efficiency (%) for 7r^Zc(3885) (Zc(3885) —>■ D^D* ), D* — > 

TT ~—>■ i, Z5° —>■ j, where i and j are described in the caption of Table imi 


s ■ x r>Tn m 1 • o Kinematic fit MC statistics Cross feed Total 

{i, BID fracking vr 4 23 GeV 4.26 GeV 4.23 GeV 4.26 GeV 4.23 GeV 4.26 GeV 4.23 GeV 4.26 GeV 


{ 0 , 0 } 

4 

5 

0 

0.6 

0.5 

0.6 

0.6 

0.2 

0.2 

6.5 

6.5 

{ 0 , 1 } 

4 

5 

3 

0.6 

0.3 

0.8 

0.8 

0.1 

0.1 

7.1 

7.1 

{ 0 , 2 } 

6 

7 

0 

0.7 

1.2 

0.8 

0.8 

0.1 

0.3 

9.3 

9.3 

{0, 3} 

6 

7 

3 

1.2 

0.9 

1.2 

1.2 

0.2 

0.0 

9.8 

9.8 

{ 1 , 0 } 

4 

5 

3 

0.5 

0.6 

0.8 

0.8 

0.1 

0.2 

7.1 

7.1 

( 1 , 1 } 

4 

5 

6 

0.7 

0.5 

1.2 

1.2 

0.1 

0.0 

8.9 

8.9 

( 1 , 2 } 

6 

7 

3 

0.9 

0.4 

1.2 

1.2 

0.2 

0.1 

9.8 

9.8 

(1, 3} 

6 

7 

6 

0.8 

0.6 

2.1 

2.1 

0.1 

0.0 

11.2 

11.2 

( 2 , 0 } 

6 

7 

0 

0.7 

0.8 

0.8 

0.8 

0.2 

0.1 

9.3 

9.3 

( 2 , 1 } 

6 

7 

3 

0.6 

0.5 

1.1 

1.1 

0.1 

0.1 

9.8 

9.8 

( 2 , 2 } 

8 

9 

0 

1.3 

1.1 

1.1 

1.1 

0.0 

0.0 

12.2 

12.1 

(2, 3} 

8 

9 

3 

0.5 

1.1 

2.0 

2.1 

2.0 

2.9 

12.7 

13.0 

{3, 0} 

6 

7 

3 

0.8 

0.6 

1.2 

1.2 

0.1 

0.3 

9.8 

9.8 

(3, 1} 

6 

7 

6 

0.6 

0.9 

2.0 

2.1 

0.0 

0.1 

11.2 

11.2 

(3, 2 } 

8 

9 

3 

1.0 

1.6 

2.1 

2.1 

2.4 

2.5 

12.8 

12.9 

(3, 3} 

8 

9 

6 

0.9 

1.0 

5.4 

5.8 

0.0 

0.0 

14.5 

14.7 
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TABLE VII. The systematic uncertainties for signal efficiency (%) for 7r+Zc(3885) (Zc(3885) —>■ D D*^), D*° 

•n^D^, D~ —>■ i, —>■ j, where i and j are described in the caption of Table HVl 


/ • -r DTri m 1 ■ 0 li'O Kinematic fit MC statistics Cross feed Total 

{i,]\ BID fracking TT Kg ^ ^3 4.26 GeV 4.23 GeV 4.26 GeV 4.23 GeV 4.26 GeV 4.23 GeV 4.26 GeV 


{A, 0} 

5 

6 

0 

0 

0.3 

0.4 

0.6 

0.6 

0.3 

0.4 

7.9 

7.9 

{B, 0} 

5 

6 

3 

0 

0.1 

0.1 

1.0 

1.0 

0.3 

0.2 

8.4 

8.4 

{C, 0} 

3 

4 

0 

4 

0.2 

0.2 

0.6 

0.6 

0.4 

0.3 

6.5 

6.4 

{D, 0} 

3 

4 

3 

4 

0.4 

0.3 

0.9 

0.9 

0.2 

0.2 

7.1 

7.1 

{E, 0} 

5 

6 

0 

4 

0.7 

0.5 

0.8 

0.8 

0.1 

0.1 

8.8 

8.8 

{F, 0} 

5 

6 

0 

0 

0.4 

0.3 

0.7 

0.7 

0.5 

0.5 

7.9 

7.9 

1} 

5 

6 

3 

0 

0.3 

0.6 

0.9 

0.9 

0.1 

0.1 

8.4 

8.4 

{B, 1} 

5 

6 

6 

0 

0.3 

0.7 

1.5 

1.5 

0.1 

0.1 

10.0 

10.0 

{C, 1} 

3 

4 

3 

4 

0.4 

0.3 

0.9 

0.9 

0.2 

0.2 

7.1 

7.1 

{D, 1} 

3 

4 

6 

4 

0.2 

0.1 

1.4 

1.4 

0.2 

0.1 

8.9 

8.9 

{E, 1} 

5 

6 

3 

4 

0.9 

0.8 

1.3 

1.3 

0.3 

0.5 

9.4 

9.4 

{F, 1} 

5 

6 

3 

0 

0.6 

0.4 

1.1 

1.0 

0.2 

0.3 

8.5 

8.4 

{A 2} 

7 

8 

0 

0 

0.6 

1.0 

0.9 

0.9 

0.2 

0.1 

10.7 

10.7 

{B, 2} 

7 

8 

3 

0 

0.5 

0.4 

1.4 

1.4 

0.1 

0.3 

11.1 

11.1 

{C, 2} 

5 

6 

0 

4 

0.4 

0.6 

0.8 

0.8 

0.2 

0.0 

8.8 

8.8 

{D, 2} 

5 

6 

3 

4 

0.4 

0.3 

1.3 

1.3 

0.2 

0.2 

9.4 

9.4 

{E, 2} 

7 

8 

0 

4 

1.0 

1.2 

1.3 

1.3 

0.0 

0.0 

11.5 

11.5 

{F, 2} 

7 

8 

0 

0 

1.0 

0.9 

1.0 

1.0 

0.3 

0.3 

10.7 

10.7 

{A, 3} 

7 

8 

3 

0 

0.9 

1.0 

1.4 

1.5 

1.2 

2.0 

11.2 

11.4 

{B, 3} 

7 

8 

6 

0 

0.0 

0.3 

2.9 

3.0 

0.0 

0.0 

12.5 

12.6 

{C, 3} 

5 

6 

3 

4 

1.1 

0.7 

1.4 

1.4 

0.3 

1.0 

9.4 

9.5 

{D, 3} 

5 

6 

6 

4 

0.8 

1.2 

2.5 

2.6 

0.0 

0.0 

10.9 

11.0 

{E, 3} 

7 

8 

3 

4 

1.2 

1.6 

2.4 

2.5 

0.0 

0.0 

12.1 

12.1 

{F, 3} 

7 

8 

3 

0 

0.8 

0.9 

1.7 

1.7 

0.2 

0.2 

11.2 

11.2 


TABLE VIII. The systematic uncertainties for signal efficiency (%) for 7r"'"Zc(3885) (Zc(3885) —f D°D* ), D* —>■ 

'k°D~ , D~ —>■ i, —>■ j, where i and j are described in the caption of Table HVl 


r- 1 DTT^ m 1 ■ „o t.'O Kinematic fit MC statistics Cross feed Total 

{i,ji BiU fracking TT Kg ^ ^3 4 25 GeV 4.23 GeV 4.26 GeV 4.23 GeV 4.26 GeV 4.23 GeV 4.26 GeV 


{A, 0} 

5 

6 

0 

0 

0.7 

0.5 

0.7 

0.7 

0.5 

0.5 

7.9 

7.9 

{B, 0} 

5 

6 

3 

0 

0.4 

0.2 

1.0 

1.0 

0.2 

0.3 

8.4 

8.4 

{C, 0} 

3 

4 

0 

4 

0.2 

0.3 

0.6 

0.6 

0.2 

0.3 

6.4 

6.4 

{D, 0} 

3 

4 

3 

4 

0.2 

0.2 

0.9 

0.9 

0.3 

0.3 

7.1 

7.1 

{E, 0} 

5 

6 

0 

4 

1.0 

0.8 

0.8 

0.8 

0.3 

0.3 

8.9 

8.9 

{F, 0} 

5 

6 

0 

0 

0.4 

0.5 

0.7 

0.7 

0.4 

0.5 

7.9 

7.9 

1} 

5 

6 

3 

0 

0.8 

0.6 

0.9 

0.9 

0.2 

0.2 

8.5 

8.4 

{B, 1} 

5 

6 

6 

0 

0.5 

0.3 

1.5 

1.5 

0.1 

0.0 

10.0 

10.0 

{C, 1} 

3 

4 

3 

4 

0.4 

0.5 

0.9 

0.9 

0.1 

0.2 

7.1 

7.2 

{D, 1} 

3 

4 

6 

4 

0.4 

0.2 

1.5 

1.4 

0.1 

0.1 

8.9 

8.9 

{E, 1} 

5 

6 

3 

4 

0.8 

0.9 

1.3 

1.3 

0.6 

0.4 

9.4 

9.4 

{F, 1} 

5 

6 

3 

0 

0.6 

0.7 

1.1 

1.0 

0.3 

0.3 

8.5 

8.5 

{A 2} 

7 

8 

0 

0 

0.8 

0.9 

0.9 

0.9 

0.2 

0.2 

10.7 

10.7 

{B, 2} 

7 

8 

3 

0 

1.1 

0.5 

1.4 

1.5 

0.2 

0.2 

11.2 

11.2 

{C, 2} 

5 

6 

0 

4 

0.8 

0.8 

0.8 

0.8 

0.2 

0.1 

8.9 

8.9 

{D, 2} 

5 

6 

3 

4 

0.6 

0.4 

1.3 

1.4 

0.3 

0.3 

9.4 

9.4 

{E, 2} 

7 

8 

0 

4 

1.4 

1.2 

1.3 

1.3 

0.0 

0.0 

11.5 

11.5 

{F, 2} 

7 

8 

0 

0 

1.1 

1.1 

1.0 

1.0 

0.3 

0.3 

10.7 

10.7 

{A, 3} 

7 

8 

3 

0 

1.1 

1.1 

1.5 

1.5 

1.4 

1.8 

11.3 

11.3 

{B, 3} 

7 

8 

6 

0 

1.3 

0.1 

2.9 

2.9 

0.0 

0.0 

12.6 

12.6 

{C, 3} 

5 

6 

3 

4 

0.8 

0.8 

1.4 

1.4 

0.2 

0.9 

9.4 

9.5 

{D, 3} 

5 

6 

6 

4 

0.1 

0.4 

2.5 

2.6 

0.0 

0.0 

10.9 

11.0 

{E, 3} 

7 

8 

3 

4 

1.6 

1.2 

2.5 

2.6 

0.0 

0.2 

12.1 

12.1 

{F, 3} 

7 

8 

3 

0 

0.6 

1.0 

1.7 

1.7 

0.2 

0.2 

11.2 

11.2 














